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ABSTRACT 


An  investigation  is  made  into  the  feasibility  of  using 
mixing-ratio  predictions  employing  multiple  linear  regres- 
sion techniques  applied  to  saturated  mixing  ratio  profiles 
at  specified  levels.   Three  such  sets  of  coefficients  are 
determined  from  Ocean  Station  Vessel  sounding  data  recorded 
for  the  months  of  January  (1967-70,  inclusive),  and  tested 
for  effectiveness  to  determine  a  "best"  prediction  scheme. 
The  resultant  "best"  set  of  coefficients  is  then  compared  to 
a  previously  designed  regression-prediction  scheme  due  to 
Weinreb  and  Crosby  (1973)  based  on  a  selection  of  island  and 
coastal  data  divided  equally  in  time  and  space.   Tests  were 
performed  to  determine  possible  sources  of  error  affecting 
the  water  vapor  specifications  made  by  the  multiple  regres- 
sion procedure. 
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I.   INTRODUCTION 

Weinreb  and  Crosby  (1973),  (henceforth  abbreviated  WC), 
presented  a  useful  technique  for  estimating  atmospheric 
moisture  profiles  to  assist  in  obtaining  more  accurate  tem- 
perature retrievals  from  satellite  soundings.   Their  estima- 
tion procedure  was  designed  to  derive  the  mixing  ratios  W  at 
the  lowest  40  levels  of  the  100-level  atmosphere  (McMillin, 
et  al,  1973)  defined  in  terms  of  linear  j-increments  of  p2'7 
given  by 

P.  =  .01(l.+(.26087836)(j-l)7/2  (1) 

The  saturated  mixing  ratio  values  W  (k)  at  11  key  levels 
are  chosen  as  predictors  from  among  the  lowest  40  j-levels 
of  eq.  (1).   In  (1),  as  j  =  01,..., 100,  P.  ranges  over 
.01,..., 1000  mb  in  the  manner  shown  in  Table  1.   Table  1 
also  indicates  the  11  fixed  k-levels  designated  for  computa- 
tional use  of  the  W  (k)-values  employed  as  predictors.   The 
data  set  used  by  WC  to  establish  the  regression  coefficients 
was  deduced  from  1100  island  and  coastal  radiosondes  with 
humidities  reported  to  200  mb.   The  WC  radiosonde  sample  was 
distributed  in  an  approximately  random  manner  with  respect 
to  season  and  latitude. 

The  forty  regression  equations  of  WC  were  presented  in 

the  form 

11 
Wj  =  W^+jj^  Cjk(Ws(k)-Ws(k));  j=l,...,40.  (2) 
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Table  1.   The  100  level  atmosphere.   Each  j  is  an  increment  of  P2'7, 
j  =  1,...100.   Asterisked  values  are  used  for  predictor 
levels  in  the  regression  analyses  conducted. 


Level 

Pressure 

Level 

Pressure 

Level 

Pressure 

j 

P(j) 

j 

P(j) 

j 

P(j) 

(mb) 

(mb) 

(mb) 

1 

.0100 

35 

30.2057 

68 

271.2453 

2 

.0225 

36 

33.0936 

69 

284.8862 

3 

.0434 

37 

36.1735 

70* 

299.0103 

4 

.0756 

38 

39.4530 

71 

313.6275 

5 

.1219 

39 

42.9395 

72 

328.7482 

6 

.1857 

40 

46.6407 

73 

344.3825 

7 

.2703 

41 

50.5643 

74 

360.5408 

8 

.3794 

42 

54.7182 

75 

377.2336 

9 

.5168 

43 

59.1104 

76 

394.4712 

10 

.6866 

44 

63.7487 

77* 

412.2642 

11 

.8928 

45 

68.6414 

78 

430.6232 

12 

1.1398 

46 

73.7966 

79 

449.5589 

13 

1.4322 

47 

79.2225 

80 

469.0820 

14 

1.7746 

48 

84.9277 

81 

489.2033 

15 

2.1718 

49 

90.9204 

82* 

509.9337 

16 

2.6288 

50 

97.2092 

83 

531.2841 

17 

3.1507 

51 

103.8027 

84 

553.2655 

18 

3.7426 

52 

110.7097 

85 

575.8889 

19 

4.4100 

53 

117.9389 

86 

599.1655 

20 

5.1584 

54 

125.4991 

87 

623.1065 

21 

5.9933 

55 

133.3993 

88 

647.7231 

22 

6.9206 

56 

141.6484 

89 

673.0266 

23 

7.9461 

57 

150.2557 

90* 

699.0285 

24 

9.0758 

58 

159.2302 

91 

725.7401 

25 

10.3158 

59 

168.5812 

92* 

753.1729 

26 

11.6723 

60 

178.3181 

93 

781.3385 

27 

13.1517 

61 

188.4501 

94* 

810.2485 

28 

14.7804 

62 

198.9868 

95 

839.9147 

29 

16.5049 

63* 

209.9378 

96* 

870.3486 

30 

18.3920 

64 

221.3126 

97 

901.5623 

31 

20.4283 

65 

233.1209 

98* 

933.5674 

32 

22.6209 

66* 

245.3726 

99 

966.3760 

33 

24.9765 

67 

258.0774 

100* 

1000.0000 

34 

27.5024 
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Equation  (2)  gives  the  format  of  the  prediction  equation 

where  j  is  the  level  of  the  mixing  ratio  prediction,  k  is 

the  predictor  level  and  W.,  W  (k)  the  respective  sample-means 

J    s 

of  W.,  W  (k)  from  the  1100  soundings.   As  seen  in  Table  1, 
J    s 

the  arrangement  of  the  predictor  levels  (denoted  by  aster- 
isks) takes  primary  advantage  of  the  specification  of  the 
W  -profile  in  the  lower  levels,  yet  allows  for  additional 
upper  level  influences. 

The  estimation  method  presented  by  WC  demonstrated  a 
strong  multiple  correlation  between  the  predicted  W.  values 
and  the  indicated  saturated  mixing  ratio  as  predictands. 
However,  in  the  collection  of  1100  sample  soundings,  strati- 
fication of  the  regression  results  of  WC  by  season  was  not 
considered.   Additionally,  in  application  to  the  computation 
of  precipitable  water  vapor  for  radiometric  (VTPR)  calcula- 
tions, island  and  coastal  radiosondes  were  considered  as 
being  representative  of  open  ocean  areas. 

In  this  thesis  several  multiple  regression  estimation 
formulas  for  water  vapor  porfiles  are  determined  based  upon 
the  W.  and  W  (k)  data  samples  from  only  ten  Ocean  Station 
Vessels  (OSV)  distributed  over  the  Atlantic  and  Pacific. 
However,  the  OSV  data  were  based  only  on  soundings  from  the 
month  of  January  in  the  years  1967  through  1970,  and  so 
should  possibly  be  considered  representative  of  midwinter 
moisture  predictions  over  the  oceans.   The  regression  formu- 
las derived  were  based  upon  1081  dependent  data  samples  and 
then  applied  for  independent  test  purposes  upon  a  smaller 
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data  sample  drawn  from  the  same  set  of  OSV-soundings .   Eval- 
uation of  the  various  estimation  methods  derived  here  was 
conducted  in  order  to  arrive  at  a  "best"  set  of  mixing-ratio 
specification  equations.   This  prediction  scheme  was  then 
compared  with  the  WC  prediction  equations  for  the  purpose  of 
revealing  the  relative  specification  significance  in  the 
equations  derived  here.   This  evaluation  was  conducted  using 
the  smaller  independent  OSV  data  sample.   Finally,  a  data- 
sample  developed  from  160  island  and  coastal  soundings,  col- 
lected from  1-28  March  1973  and  similar  to  those  used  by  WC 
(1973),  was  subjected  to  the  "best"  specification  system  de- 
veloped here  and  compared  with  independent  testing  based  on 
the  use  of  the  WC  specification  system. 

A  considerable  amount  of  linear  interpolation  is  required 
to  adjust  sounding  data  to  the  100-level  atmosphere,  specif- 
ically to  the  levels  61  through  100.   In  this  thesis  and  in 

the  work  done  by  WC  simple  linear  interpolation  on  pressure 

t  hs 
raised  to  the  2/7    power  as  listed  in  Table  1  is  used.   It 

is  reasonable  to  assume  that  the  regression  coefficients  de- 
termined by  WC  are  not  independent  of  one  another  over  the 
40-level  system,  and  so  could  be  reduced  to  a  reasonably 
concise  but  smaller  predictand  set,  say  J=l,...,19.   For 
purposes  of  this  thesis,  the  best-fit  multiple  regression 
equation  in  the  form  of  (2)  was  therefore  applied  only  at 
the  19  levels  listed  in  Table  2.   It  is  to  be  noted  that 
this  set  of  19  levels  is  a  subset  of  the  40-level  predictand 
set  of  WC  and  that  it  also  includes  the  11  levels  used  as 
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Table  2.   19-level  atmosphere.   The  pressure  levels  P(J)  are 
obtained  by  linear  scaling  of  P2'7over  the  interval 

188.45  mb  to  1000  mb  in  19  increments.   Column  3 

indicates  the  subset  of  predictor  levels. 

Predictand       W  (k)  Predictor 

s 

Level  Pressure  Level         Levels 

(J)  P(J)  (mb)            (k) 

1  188.45 

2  209.94  1 

3  245.37  2 

4  271.24 

5  299.01  3 

6  328.75 

7  360.54 

8  412.26  4 

9  469.08 

10  509.93  5 

11  553.26 

12  .     599.16 

13  647.72 

14  699.03  6 

15  753.17  7 

16  810.25  8 

17  870.35  9 

18  933.56  10 

19  1000.00  11 
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predictor-levels  in  their  system.   It  was  assumed  that  re- 
stricting Eq.  (2)  to  19-predictand  levels  would  result  in 
little,  if  any,  loss  of  resolution  of  the  estimated  mixing 
ratio  profiles  over  the  lowest  40  WC  predictand  j-levels. 
Hence,  comparison  of  the  regression  results  deduced  here 
will  be  with  those  results  of  WC  at  the  19-predictand  levels 
adopted  for  use  in  this  study. 

Regression  analysis  was  carried  out  on  the  IBM  360  com- 
puter using  the  BMD03R  statistical  program  (Dixon,  1973). 
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II.   DATA  PROCESSING 

A.  OCEAN  STATION  VESSEL  RADIOSONDES 

The  raw  data  used  in  this  study  consisted  of  a  set  of 
January  radiosondes  containing  temperature  and  relative  hu- 
midity data  in  accordance  with  standard  radiosonde  reporting 
procedures,  i.e.  available  mandatory  level  plus  any  signifi- 
cant levels  (Radiosonde  Observations  Handbook,  1969),  from 
the  10  OSVs.   Table  3  gives  the  approximate  latitude  and 
longitude  for  each  of  the  OSVs  in  addition  to  the  sample- 
size  contributed  by  each  to  the  total  data  set.   Also  listed 
in  Table  3  is  a  similar  breakdown  of  the  sample  size  contri- 
butions to  the  total  dependent  and  independent  data  samples 
(abbreviated  as  OSV-1,  and  OSV-2,  respectively).   The  latter 
two  data  samples  will  be  discussed  later  in  this  section. 

Each  sounding  included  in  the  OSV  data  samples  was  re- 
quired to  have  temperature  data  up  to  at  least  150  mb  and 
relative  humidity  data  to  450  mb.   Both  mandatory  and  sig- 
nificant levels  were  combined  to  produce  temperature  and 
relative  humidity  profiles  at  pressure  levels  arranged  in 
order  from  1000  mb  to  150  mb. 

B.  REDUCTION  OF  RADIOSONDE  OBSERVATIONS 

Once  a  radiosonde  was  arranged  in  a  pressure-decreasing 
order  the  Goff-Gratch  formula  in  the  form  (List,  1963) 
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Table  3.   Radiosonde-data  sample  sizes  from  Ocean  Station  Vessels,  and 
their  geographic  locations. 


POSITION 
OSV   LATITUDE  LONGITUDE 


NUMBER  OF  SOUNDINGS  IN  SAMPLES 
TOTAL   DEPENDENT  DATA   INDEPENDENT  DATA 
OSV       (OSV-1)  (OSV-2) 


A 

62.  5N 

33. 0W 

194 

97 

B 

56.  5N 

50.  5W 

154 

77 

C 

53.  5N 

35.  5W 

189 

95 

I 

59.  ON 

19.  0W 

212 

106 

J 

53.  5N 

19.  5W 

206 

103 

P 

50.  ON 

145. OW 

202 

101 

D 

44.  5N 

41.  OW 

169 

113 

E 

35.  ON 

48.  5W 

185 

123 

N 

30.  ON 

140. OW 

205 

136 

V 

38.  ON 

165. OW 

194 

130 

97 

77 

94 

106 

103 

101 

56 

62 

69 

64 


TOTALS 


1910 


1081 


829 
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log(e  /Po)  =  -7.90298((T  /T)-l)+5. 02808  log(T  /T) 
s  s  s 


-1.3816  x  io-7(1011-33,t(1~T/Ts)-i) 

+8.1328  x  10-3[10"3-lt9149((Ts/T)~1)-l],     (3) 


where  T  =  373.16°K,  P0  =  1013.246  mb,  and  T  =  T(P)  is  the 
temperature  at  pressure  P,  was  used  to  determine  saturated 
vapor  pressures  W  .   From  the  value  of  e  ,  W  (P)  was  deter- 

o  S     S 

mined  from 

.622  e 
WS(P)  -  p-~  (4) 


which  in  turn  yielded  the  actual  mixing  ratio  W(P)  upon  appli- 
cation of  the  relationship  between  RH  (as  a  fraction)  and  W 

W(P)  =  RH  x  w  (P)  (5) 

Following  this  procedure,  simple  linear  interpolation  of 
both  T(P)  and  W(P)  relative  to  the  [P(J)]2/7  values  of  Table 
2  was  performed,  yielding  a  temperature  profile  T(J)  up  to 
188.45    mb  and  a  mixing  ratio  profile  W(J)  to  469.08   mb. 
Equations  (3)  and  (4)  were  again  needed  to  establish  the 
11-predictor  values  of  W  (k);  in  this  case  T  =  T(k),  P  =  P(k), 
k=l,...,ll,  as  shown  in  Table  2. 

Since  mixing  ratio  (or  RH)  profiles  were  usually  unavail- 
able in  the  radiosondes  above  450  mb,  the  well  known  cubic 
power-law  profile  of  Smith  (1966)  expressing  W  in  the  form 

W  =  W   .(P/P   ,)X,    A  =  3  (6) 

refv  '  ref 
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was  used.   Here  W   _  and  P   _  signify  the  simultaneous  con- 

ref       ref 

ditions  at  the  reference  level,  taken  here  as  P   -  =  469.08  mb 

ref 

Equation  (6)  was  employed  to  compute  "observational"  values 

of  W  at  all  levels  above  P   .to  complete  the  19  levels  of 

ref       r 

Table  2. 

Completion  of  the  above  procedures  through  (6)  yielded, 
for  each  sounding,  an  array  of  19  values  of  W  and  11  values 
of  1  .   These  values  were  arranged  in  a  final  30-variable 
format  [W( 1) , W(2 ) , . . . , W( 19) , Wg( 1 ) , Wg(2 ) , .. . ,Ws(ll)].   This 
format  was  found  to  be  most  convenient  accessing  the  least 
squares  statistical  program,  BMD03R. 

C.   DIVISION  OF  DATA  SAMPLES 

The  collection  of  1910  OSV  soundings  was  divided  into 
two  major  data  samples,  OSV-1  consisting  of  1081  total 
soundings,  and  OSV-2,  comprised  of  the  remaining  829.   OSV- 2 
data  were  obtained  by  selecting  a  random  mix  of  every  second 
to  every  third  sounding  from  the  total  set  of  soundings 
available  from  each  ship,  with  the  remainder  going  into  the 
dependent  sample  OSV-1.   Table  3  lists  the  contributions  to 
each  sample  on  a  ship  by  ship  basis.   OSV-1  was  used  for  de- 
termination of  the  dependent-test  set  of  regression  equa- 
tions, for  W(J)  specified  in  terms  of  the  predictor-variables 
W  (1),...W  (11).   OSV-2  was  used  for  independent  or  verifi- 
cation  tests  of  regression  results  developed  using  OSV-1. 

An  additional  subdivision  of  the  OSV  data  was  made  to 
determine  latitudinal  effects,  if  any,  on  the  moisture  esti- 
mation regressions.   This  subdivision  was  suggested  after 
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consideration  of  the  distributions  of  the  mean  mixing  ratios 
W(J)  for  each  ship  station  in  the  OSV-1  data  sample. 

Table  4  lists  the  mean  W-values  by  level  for  the  entire 
dependent  sample  set,  the  subdivisions  of  the  sample  accord- 
ing to  <{>>50N,  and  cf><50N,  as  well  as  by  individual  ships.   It 
seemed  plausible  to  use  50N  as  the  dividing  point  for  the 
latitudinal  zones  as  this  arrangement  provided  reasonably 
large  sample  sizes  to  stratify  the  regression  analyses  in 
addition  to  maintaining  the  apparent  stratifications  demon- 
strated by  the  mean  distributions.   Although  Ship  Station  P 
was  located  at  50N,  its  W-profile  statistics  seemed  to  be 
described  more  closely  by  the  means  of  the  northern  set  as 
its  verification  statistics  suggest  (Sec.  IV.B.l). 

Verification  of  the  estimation  equations  derived  by  the 
latitudinal  stratification  of  the  dependent  sample  set  was 
conducted  on  the  stratifications  (a)  all  <j> ,  (b)  <J)_>50N,  (c) 
<j)<50N  as  well  as  for  the  individual  ship-data  samples  re- 
served for  the  independent  OSV-2  sample.   The  results  of  the 
tests  on  the  individual  ships  are  tabulated  in  Appendix  B. 

D.   ISLAND  AND  COASTAL  RADIOSONDE  SAMPLE 

Data  collected  under  the  auspices  of  the  NATIONAL  ENVI- 
RONMENTAL SATELLITE  SERVICE  (NESS)  during  March,  1973,  were 
used  for  additional  testing  of  the  estimation  procedures. 
This  data  set  provided  a  total  of  160  typical  radiosondes 
containing  moisture  data  from  the  surface  up  to  at  least  lev- 
el J=9  of  the  pressure  scale  of  Table  2  but  with  temperatures 
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up  to  at  least  level  1.   The  format  of  the  soundings  was 
much  the  same  as  that  of  the  OSV  data  except  that  the  mois- 
ture element  was  given  in  terms  of  the  dew  point  T.  rather 

d 

than  RH.   This  difference  required  the  use  of  Eq.  (3)  with 
e   and  T  replaced  by  e  and  T  ,  respectively,  followed  by 

S  CI 

application  of  Eq.  (4)  to  obtain  the  W-values.   All  other 
procedures  in  this  section  apply  to  the  NESS  data  to  obtain 
soundings  in  the  form  [W(l) ,W(2) , . . . ,W(J) ,W  (1),W  (2),... 

to         to 

W  (11)],  where  J  is  at  least  as  high  as  level  9.   Recall  that 
W  (1),W  (2 ),..., W  (11)  depend  simply  upon  having  the  temper- 

to  to  to 

ature  profile  complete  in  the  form  of  Table  2  to  level 
P(2)  =  209.94  mb,  and  that  for  purposes  of  computing  W  (k), 

to 

the  indexing  k=l,...,ll  denotes  the  indicated  subset  of  the 
J-values  shown  in  column  3  of  Table  2. 


E.   VERTICAL  WEIGHTING 

With  the  mixing  ratios  estimated  at  most  of  the  19  levels 
of  Table  2,  it  was  possible  to  compute  precipitable  water 
vapor  for  each  sounding  as  a  measure  of  the  vertically  inte- 
grated mixing  ratios.   The  precipitable  water  vapor,  calcu- 
lated by  means  of  the  finite-difference  integration,  Eq.  (7). 

1<J<9 

*  (wj+i+wj)(pj+i-pj) 

TT    1000   J=18   J  +  1   J    J  +  1   J  m 

U     980  2 

provided  a  basis  for  testing  "estimated"  precipitable  water 
vapor  against  the  corresponding  "observed"  U-values  deduced 
from  the  W-profiles  summed  over  the  same  atmospheric  levels 
from  level  J  downward.   In  the  case  of  OSV  soundings  the 
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summation  of  Eq.  (7)  extended  up  to  J=l;  however  in  the  case 
of  the  NESS  data,  the  top  level  of  humidity  reports  extended 
up  to  a  level  J  in  the  range  1<J<9. 
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III.   STATISTICAL  CONCEPTS  AND  PROCEDURES 

A.   MULTIPLE  LINEAR  REGRESSION 

This  thesis  makes  use  of  multiple  linear  regression 
(MLR)  techniques  to  develop  an  effective  specification 
scheme  to  estimate  mixing  ratio  values  at  the  19  previously 
designated  pressure  levels,  using  as  predictors,  the  satu- 
rated mixing  ratios  obtained  at  the  11  key  levels  specified 
by  WC.   In  the  discussion  to  follow,  the  phraseology  "de- 
pendent variable"  refers  to  observed  mixing  ratios  to  be 
specified  in  terms  of  the  independent  variables  when  MLR 
techniques  are  to  be  applied.   The  independent  variables  are 
always  the  set  of  saturated  mixing  ratios  computed  at  the  11 
specified  levels  in  the  vertical  soundings. 

Assuming  that  the  dependent  variable  Y  depends  on  k  in- 
dependent variables  (X- ,X  , . . . , X.) ,  and  that  there  is  a  sam- 
ple of  N  observations  of  the  k+1  variables  (Y, X- ,X  , . . . ,Xfc) , 
MLR  techniques  can  be  applied  to  find  the  best-fit  multiple 
regression  coefficients  C..   Y  can  be  estimated  by  the  re- 
gression equation 

k 
Y  =  A  +   £   C .X .  (8) 

where  the  superscripted  symbol  Y  denotes  "estimation  of  Y." 

Also  it  follows,  from  the  requirement  that  the  regres- 
sion plane  pass  through  the  sample  mean,  that 

k 

A  =  Y-  .L  C.X..  (9) 

J=l   J  J 
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The  quantities  (  )  signify  the  respective  means  based  on  the 
sample  size  N. 

It  should  be  noted  that  Eqs.  8  and  9  can  be  combined  to 
produce  the  form 


Y  -  Y  =  .2   C.(X.-X.) 


(10) 


where  the  C.'s  are  the  best-fit  regression  coefficients  de- 
termined  by  MLR  techniques. 

The  regression  coefficients  of  (10)  are  determined  by 
solving  the  system  of  normal  equations  given  in  determinant 
form  by 


Ci jC2  » ' * '  'Ck 


a    a 

11      12 


'kl 


a 


a.  ,   a-„  „  ...  a 

2  1      2  2 


ik 

a 

iy 

2k 

= 

a, 

2y 

kk 

a, 

ky 

^j-jl; 


where  a. .  and  a.   are  given  by 
ij       iy 


N 

a.  .  =   £  (X.  -X.  )(X.  -X.), 
ij    m=lv  lm  \'  v  jra   j 


(12) 


and 


N 


..   =   £  (X.  -X.  )(Y  -Y) 
iy    m=lv  lm   i/v  m 


(13) 


respectively . 

In  Eqs.  (12)  and  (13)  it  is  understood  that  (Y^^X^  , 

. ...xfm^)  constitute  the  mth  sample  of  Y  on  X    For  simpli- 
k  x 

fication  the  superscript  in  this  sample  designation  has  been 
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written  as  a  subscript  in  (12)  and  (13).   From  Eq.  (12), 
a. .  =  a.. ,  and  therefore,  the  matrix  of  {a. .}  in  (11)  is  sym- 
metric.  Equation  (11)  represents  a  system  of  k  equations  in 
the  k  unknowns  C.  which  can  be  solved  for  a  unique  set  of 
regression  coefficients,  provided  the  k-by-k  determinant  of 
the  left  hand  side  of  (11)  is  non-singular.   If  the  multipli- 
cation of  the  left  hand  terms  of  (11)  is  carried  out,  one 
can  obtain  k  equations  in  the  form 

ak.C,  +  ak2C2  +  •••  +  akkCk  ■  aky  (14> 

which  can  be  solved  (Crow,  et  al.,  1955)  by  matrix  inversion 
methods  applied  to  {a. .}. 

B.   SIGNIFICANCE  TESTS 

1.   On  Dependent  Regression  Equations 

The  coefficient  of  determination  of  an  MLR  technique 
is  definable  as  the  ratio  of  the  sum  of  squares  due  to  the 
regression  to  the  total  sum  of  squares.   This  parameter  may 
be  shown  to  be  the  square  of  the  multiple  correlation  coeffi- 
cient, R2, 

M 

N 
C  y,    (Y   -Y ) 2 1 
R2  _  Sum  Sq.  Explained  by  Reg.   _   i=lv  i  (15) 

M  "  Sum  Sq.  of  Total  Deviation     N   A   _ 

[i=l(Yi"Y)2] 

and  is  used  as  a  measure  of  the  efficiency  of  the  regression 
performed  by  the  BMD03R  program.   The  total  sum  of  squares 

N  N 

S.S.  total  =  ^(Y.-Y)2  =  iJ1[(Yi-Y)  +  (Y.-Y.)]2  =  (N-l)S2       (16) 
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defines  the  variance  S2    S.S.  total  may  be  expanded 

N  N  N 

S.S.  total  =  i£1(Yi-Y)2  =  iS1(Yi-Y)2+i21(Y.-Yi)2     (17) 

(a)     +    (b) 

which  is  a  combination  of  the  S.S.  due  to  regression  (a), 
and  the  residual  S.S.  (b). 

If  the  regression  were  perfect,  term  (b)  of  (17) 
would  vanish  and  R2  would  be  1;  conversely,  if  Y  were  totally 
independent  of  X.,  R2  would  be  0.   Combination  of  Eqs.  (15) 
and  (17)  shows  that  the  residual  sum  of  squares  after  appli- 
cation of  the  regression  (8)  is 

N 

Res.  S.S.  =  (1  -  R2)   £  (Y.-Y)2  . 
v      M'  i=lv  l  ' 

A  method  of   determining   the   significance   of    a  predic- 
tion  equation   having   k  predictors    is   the   F-statistic  with  k 
and  N-k-1   degrees   of    freedom,    defined  by 

_  Mean   Square  Explained  by   the  Prediction  (  *\R\ 

'  Mean   Square  Unexpl.    after   the  Prediction    ' 

or 

RM/k 
F(k, N-k-1)    =  - .  (19) 

(l-R2)/(N-k-l) 

Equation  19  is  a  direct  application  of  the  defini- 
tion (18)  of  the  F-statistic  for  the  appropriate  degrees  of 
freedom.   For  the  numerator  of  (19),  k  is  the  number  of  in- 
dependent variables  used  in  specifying  the  mean  square  ex- 
plained of  (18).   In  the  denominator  of  (19)  the  degrees  of 


29 


freedom  (df)  is  equal  to  N-k-1,  corresponding  to  the  total 
sample  size,  reduced  by  the  11  known  coefficients  and  the 
value  of  the  sample  mean.   In  the  tests  on  the  dependent 
sample  k=ll  throughout,  although  N  varies  from  test-sample 
to  sample. 

Critical  values  of  the  F-statistic,  designated  F  , 

I  D  c 

can  be  found  in  standard  texts  on  statistics  for  the  appro- 
priate degrees  of  freedom  at  specified  levels  of  signifi- 
cance.  If  the  F-statistic  of  (18)  exceeds  F, 0 i (k, N-k-1) , 
there  is  less  than  1%  probability  that  the  sample  came  from 
a  population  in  which  there  is  no  linear  relationship  be- 
tween Y  and  X..   Such  a  statistical  outcome  for  the  regres- 
sion is  described  as  corresponding  to  a  99%  level  of 
confidence. 

A  third  related  statistic  which  is  listed  as  output 
by  BMD03R  is  the  standard  error  of  estimate  (after  applica- 
tion of  the  regression)  is  defined  by 

(STD.  ERR.,2E  s-^°rs  =  z{lCi]2  <2o> 

Equation  (20)  shows  that  the  standard  error  is  simply  the 
square  root  of  the  mean  of  the  errors  squared.   The  standard 

error  is  related  to  R2  by 

M 

(N-1)S2(1-R2) 
(STD.  ERR.)2  =  (NXl)  (21) 
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The  three  statistics  defined  in  (15),  (18),  and  (20) 
were  considered  as  comparative  measures  of  the  MLR  techniques 
of  this  work  and  were  conveniently  provided  as  a  standard 
output  from  the  BMD03R  program.   Values  of  R„,  STD.  ERR., 
and  F  are  presented  in  tabular  form  in  Section  IV,  based 
upon  MLR  techniques  applied  to  relevant  subdivisions  of  both 
dependent  and  independent  samples. 

2.   On  Independent  Data  for  Level  by  Level  Analysis 

Similar  test  procedures  were  applied  to  evaluate  es- 
timates of  the  regression  methods  from  Sec.  III.A.l  to  the 
independent  OSV-2  sample.   Slight  variations  in  the  degrees 
of  freedom  were  required  in  specifying  the  efficiency  sta- 
tistics defined  by  Eqs.  (15),  (18),  and  (20)  to  the  inde- 
pendent samples.   In  the  regression  analyses  applied  to  the 
dependent  sample,  11  predictors  and  the  mean  were  used  in 
arriving  at  N-12  degrees  of  freedom  in  the  denominator. 
However  once  the  coefficients  were  applied  to  the  OSV-2  sam- 
ple, there  was  essentially  only  one  way  to  form  the  estima- 
tions Y,  using  the  prescribed  set  of  coefficients.   Therefore, 
in  this  section  the  df  corresponding  to  k  was  taken  as  1. 
R*  was  also  somewhat  modified  and  was  termed  simply  R£FF> 
defined  as  the  effective  coefficient  of  determination.   The 
Rp-p-p  values  were  defined  by 

R2    =  -,  _  S.S.  Errors  (22) 

nEFF        S.S.  Total 
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The  F-statistic  and  the  standard  error  of  estimate 
were  calculated  in  a  manner  similar  to  Eqs.  (19)  and  (20) 
with  k=l,  but  with  N-k-1  replaced  by  N-l.   The  resulting  ex- 
pression for  these  statistics 

R2 
EFF 
F  =  — (23) 

(1-R2FF)/(N-1) 

and 

E(Y  -Y. )2 
(STD.  ERR.)2  =  — N_-.   —  (24) 

were  used  when  conducting  significance  tests  on  the  inde- 
pendent data  set. 

The  statistics  were  collected  and  tabulated  at  each 
of  the  19  prediction  levels.   Comparison  of  their  values 
from  the  various  estimation  methods  provided  a  means  of  de- 
termining whether  a  latitudinally  dependent  specification 
scheme  over  the  OSV  sounding  samples  was  more  effective  than 
a  composite-latitude  specification  scheme.   After  determina- 
tion of  a  "best"  scheme,  specification  usefulness  was  com- 
pared in  the  same  manner  with  the  WC  prediction  scheme. 
(Here  WC  denotes  the  Weinreb-Crosby  regression  system. ) 
3.   Independent  Data  by  Precipitable  Water  Vapor 

Precipitable  water  vapor  values,  both  actual  and 
predicted,  calculated  by  Eq.  (7)  were  subjected  to  a  simple 
linear  regression  of  U  on  U  by  BMD03R,  in  order  to  obtain  a 
vertically  weighted  comparison  of  the  various  specifications 
The  three  statistics  described  in  Section  III.B.l  were  used 
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to  judge  the  effectiveness  of  the  estimates  by  the  various 
methods . 
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IV.   STATISTICAL  RESULTS   . 

A.   DEPENDENT  SAMPLE 

This  phase  of  the  research  consisted  of  collecting, 
tabulating,  and  evaluating  the  results  of  three  MLR  tests 
carried  out  by  BMD03R  on  the  dependent  OSV-1  sample  of  sound- 
ings.  The  first  MLR  provided  the  12  coefficients  and  the 
statistics  for  each  of  the  19  prediction  equations  to  be 
applied  over  the  composite  sample  of  ship-soundings.   The 
second  and  third  MLR's  provided  similar  sets  of  coefficients 
and  statistics  for  the  latitudinal  divisions  of  the  soundings 
for  (J)>_50N  and  <|><50N,  respectively.   These  three  sets  of  re- 
gression equations,  employed  in  the  form  of  Eq.  (8),  will  be 
referred  to  as: 

1)  COEF-1  for  the  one  over-all  set,  regardless  of  cf> , 

2)  C0EF-2N,S  using 

a)  The  coefficient  set  C0EF-2N  when  <J>>_50N,  and 

b)  C0EF-2S  when  cj><50N, 

3)  The  composite  set  of  COEF-2  without  the  N  or  S 
designation. 

The  latter  set  implies  a  composite  estimation  by  one  or 

other  of  the  two  components  of  the  applicable  equations  in 

the  appropriate  latitudinal  zones.   The  sample  size  is  the 

same  as  for  COEF-1. 

Table  5  lists  the  values  of  Rw,  the  standard  error  of 

M 

estimate,  and  F  as  obtained  from  BMD03R  for  each  MLR  experi- 
ment at  each  of  the  19  levels  listed  in  Table  2.   The  most 

34 


o 

o 

o 

o 

o 

o 

o 

O 

o 

o 

O 

o 

o 

O 

o 

o 

o 

o 

o 

rs 

>fr 

rn 

ro 

o 

ro 

rH 

o 

-^ 

>t 

■—4 

in 

r- 

ro 

O 

s0 

CO 

co 

sf 

o 

CNI 

r~v 

fa 

1 

CM 

r\j 

CNJ 

CNJ 

CNJ 

ro 

a 

vt 

LO 

«— ( 

^O 

r-4 

00 

O 

CO 

CNJ 

>$- 

r«- 

O 

vl- 

Cn 

O 

CNJ 

CNJ 

CNJ 

CNJ 

CNJ 

CNJ 

ro 

ro 

CO 

CNJ 

CNJ 

r-t 

CNJ 

CNJ 

ro 

>t 

o 

lf\ 

-o 

o 

UO 

r-4 

CNI 

o 

II 

* 

>_• 

iH 

>— 1 

o 

O 

0> 

CNJ 

in 

O 

f— 4 

o 

o 

O 

•J- 

o 

CNJ 

ro 

O 

o 

ro 

o 

1 

w 

o 

CO 

•=- 

CO 

ex) 

CO 

CNJ 

O 

CO 

r~- 

<f 

ro 

CO 

o 

ro 

o 

r-4 

O 

X 

o 

CO 

in 

CM 

K 

m 

rs. 

p- 

p- 

CO 

CO 

CNJ 

m 

o 

v*- 

o 

r- 

>*- 

vO 

CNJ 

in 

o 

K 

h- 

CNJ 

1 

in 

in 

in 

in 

m 

o 

vO 

nO 

O 

-£> 

m 

in 

in 

sO 

vO 

fs. 

!*- 

CO 

o 

o 

fa 

w 

0) 
4J 
CD 
CO 

o 

CJ 

<NJ 

ro 

ro 

CNJ 

co 

CNJ 

CO 

•—4 

00 

O 

in 

vO 

|S» 

CO 

ro 

-G 

00 

CO 

o 

T5 

u 

in 

rH 

ro 

in 

co 

sJ- 

V— 1 

r- 

CNJ 

vO 

CNJ 

in 

o 

ON 

in 

ir, 

ON 

ro 

■J- 

-P 

u 

4-1 

CO 

fa 

-a 

CNJ 

ro 

^ 

in 

00 

CNJ 

<NJ 

N- 

ON 

r- 

CNJ 

o 

CNJ 

sO 

in 

£> 

«;- 

0> 

c 

CD 

•—I 

Cv) 

ro 

>t 

in 

00 

o 

—4 

CNJ 

CNI 

t— 4 

O 

ON 

•H 

CJ 

•H 
U-l 

IW 

0) 

• 

O 

o 

CO 
CD 

o 

O 

o 

o 

o 

o 

o 

o 

O 

o 

o 

o 

o 

O 

o 

o 

O 

o 

o 

iH 

«*■ 

o 

ro 

CNJ 

ro 

o 

CO 

r- 

r\l 

—* 

CNI 

f- 

ro 

in 

in 

o 

in 

CNJ 

CO 

OS 

a 

•~N 

fa 

(J 

e 
co 

CO 

CD 

CNJ 

er, 

ro 

ro 

ro 

CNJ 

<h 

CNJ 

sO 

Nt 

CNJ 

CNJ 

r- 

o 

vO 

o 

CNJ 

rr> 

r— 

^ 

t> 

co 

CO 

co 

ro 

ro 

CO 

ro 

•J" 

«*• 

-* 

vt 

ro 

ro 

ro 

««*• 

■D 

O 

0> 

ro 

u— i 

lO 

r-4 

rH 

CO 

o 

60 

II 

4-) 

•H 

c 
a) 
E 
P. 

13 

3 

CNJ 

& 

ro 

CNJ 

CNJ 

r-4 

vO 

•—4 

in 

r-4 

a« 

ro 

r- 

in 

CO 

|s- 

IS. 

CO 

o 

S 

S 

—4 

<T 

vO 

O 

vO 

<T 

«*■ 

lO 

N. 

O 

o 

ro 

is. 

00 

co 

v0 

m 

f-i 

o 

O 

CM 

fa 

CNJ 

cj 

CNJ 

CNJ 

<NJ 

CNJ 

ro 

r— 

oo 

r-- 

r— 

CNJ 

Nt 

in 

CO 

IT. 

— • 

ON 

o 

o 

CO 

| 

»o 

-o 

mO 

vJ-> 

vjO 

sO 

vO 

vO 

vfi 

vaO 

nO 

vi5 

mO 

.o 

vO 

fs- 

CO 

co 

o 

r-l 

rH 

fa 

> 

1 

fa 

CD 

> 
C/3 

o 

o 

ON 

vO 

in 

■* 

•J- 

CNJ 

O 

ON 

>*■ 

ro 

r- 

>*• 

CO 

in 

00 

m 

co 

-H 

O 

T3 

u 

CO 

r-4 

CO 

in 

ro 

o 

f-1 

in 

ro 

in 

r— 

cr- 

is. 

ro 

on 

ro 

m 

r-l 

o 

60 

c 

•H 
,4 

4J 

c 

<D 

-P 
CO 

u 
fa 

o 

—4 

-H 

CNJ 

CO 

>d- 

m0 

o> 

O 

CNJ 

o 

o 

o 

vt 

o 

CNJ 

ON 

O^ 

s* 

— i 

eg 

CNJ 

vj- 

in 

vO 

r- 

is. 

vO 

in 

sO 

3 

T3 

13 

TJ 

CD 

0) 

c 

TJ 

•H 

CO 

E 

o 

<J" 

o 

O 

O 

O 

o 

o 

o 

O 

O 

O 

o 

o 

o 

O 

o 

o 

o 

4-> 

o 

Nf 

•—4 

CO 

>o 

CNJ 

>rf- 

vt 

a^ 

in 

CNJ 

r~ 

o 

r-4 

ro 

ro 

•4- 

r-4 

<r- 

r— 

O 

14-1 

^ 

fa 

r— 1 

m 

ro 

v1" 

m 

<*■ 

r~ 

vO 

vf 

r—l 

O 

CO 

—4 

•d- 

<*■ 

r-4 

ro 

in 

o 

co 

OS 

rH 

r— 

p- 

p- 

r^ 

p* 

CO 

(T> 

o 

o-* 

CO 

r» 

m 

^0 

p- 

O 

O 

0^ 

r- 

in 

o 

co 

CO 

^^ 

r- 1 

(Nt 

m 

IS- 

•H 

O 

o 

4J 
CO 

§ 

rH 

•H 

CQ 

II 

4J 

S 

co 

CNJ 

p- 

o 

in 

r- 

«rf- 

CO 

— -i 

r- 

<7\ 

J- 

■& 

CNJ 

CO 

-H 

vO 

o 

o> 

CO 

>. 

3 

o 

in 

in 

cr- 

o 

i— ( 

p* 

ro 

CNJ 

vj- 

00 

ro 

—4 

CO 

r-4 

o 

UD 

vh 

o 

4-1 

XI 

cxT 

m 

m 

in 

in 

o 

co 

O 

CNJ 

o 

r- 

vt 

r-4 

CNJ 

IT. 

o 

o 

vO 

fNJ 

•j- 

CO 

C/3 

H 

>o 

o 

o 

-o 

vO 

sO 

r— 

p- 

r»- 

O 

O 

O 

^O 

>o 

p- 

r» 

CO 

0^ 

cr- 

14-1 

CN 

1 

o 

1 

fa 
fa 

>> 

w 

o 

. 

, 

CNJ 

r-t 

CO 

>*■ 

-*• 

o 

in 

r- 

m 

ro 

is. 

r-t 

O 

co 

4- 

0> 

co 

st 

p- 

r-l 

CO 

o 
o 

u 

-a 

u 
u 

CNJ 

p- 

Nt> 

vO 

f- 

r- 

>*• 

& 

r-4 

r— 

r- 

CO 

CO 

00 

p* 

o 

o 

0N 

CNJ 

g 

•x) 

CO 

fa 

r— 4 

.—l 

CNJ 

.-o 

^ 

vO 

ON 

O 

00 

r». 

ro 

^J" 

CO 

r-4 

r-4 

r-4 

in 

>* 

O 

3 

c 

—t 

CNJ 

rO 

>*• 

nO 

r- 

ON 

o 

O 

ON 

00 

CO 

CO 

• 

CO 

rH 

T-l 

m 

<U 

1—1 

rO 

CD 

> 

r-4 

CNJ 

ro 

>* 

in 

O 

r^ 

00 

0^ 

O 

r-4 

CNJ 

(O 

«d- 

in 

O 

r- 

co 

0> 

CO 

CD 

H 

fa 

35 


stringent  critical  F-value,  F  (k,N-k-l),  for  the  entire  dis- 
play at  the  1%  significance  level  occurs  for  the  degrees  of 
freedom  indicated  by 

F  (11,490)  =  2.98, 

which  occurs  in  the  southern  subdivision  of  OSV-1.   As  can 
be  seen  from  Table  5,  all  F-values  exceed  this  critical  val- 
ue by  a  rather  impressive  amount.   The  inference  was  drawn 
that  each  estimation  procedure  so  determined  was  based  on  a 
strong  linear  dependence  of  W  on  W  (k).   Additionally,  com- 
parison  of  the  standard  errors  of  Table  5  with  the  corre- 
sponding OSV-1  sample  standard  deviations  (listed  in  Table 
4)  indicates  a  significant  reduction  of  the  variances  by  the 
regression  estimates. 

An  example  of  a  regression  equation  at  level  J=19  (cor- 
responding to  1000  mb)  is  given  by 

W(19)  =  2. 1082+0. A979  W  (l)-0.5900  W  (2)-1.1938  W  (3) 

s  s  s 

-0.1600  W  (4)+0.2194  W  (5)-0.2095  W  (6)-0.1151  W  (7) 
s  s  s  s 

+0.2814  W  (8)+0.0815  W  (9)+0.9076  W  (10)-0.0313  W  (11)     (25) 
s  s  s  s 

Equation  (25)  is  one  prediction  equation  of  the  set  of  19 
equations  developed  from  COEF-1.   A  listing  of  the  three  co- 
efficient sets  comprising  COEF-1,  C0EF-2(N,S)  as  well  as  the 
WC  coefficient  set  is  contained  in  Appendix  A.   It  should  be 
noted  that  the  first  term  on  the  right  side  of  (25)  results 
from  Eq.  (9)  with  all  the  sample  means  known. 
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Examination  of  Table  5  does  not  reveal  convincing  evi- 
dence with  which  to  determine  a  "best"  prediction  method  be- 
tween COEF-1  and  COEF-2  as  all  F-values  demonstrate  at  least 
a  1%  significance  level  for  dependent  W-specif ications .   For 
purposes  of  this  study,  a  highly  significant  improvement  in 
specifications  by  COEF-2  over  COEF-1  would  be  required.   This 
requirement  was  deemed  essential  to  override  the  added  com- 
puter storage  and  procedures  required  to  implement  a  double 
set  of  coefficients.   Clearly,  Table  5  does  not  offer  this 
evidence. 

B.   INDEPENDENT  SHIP  SAMPLE 

1 .   An  Illustration  of  the  Usefulness  of  Lat itudinally- 
Stratified  Coefficient-Set 

Recall  that  in  Section  II,  the  question  of  a  latitu- 
dinal stratification  of  the  water  vapor  prediction  equations 
was  raised.   The  reason  for  considering  this  approach  is 
well  illustrated  by  the  level-by-level  results  for  Ship  P 
presented  in  Table  6  (based  here  upon  independent  test  data). 
The  results  of  Table  6,  e.g.  larger  values  of  REFF  resulting 
from  the  C0EF-2N  regression  scheme,  indicated  stronger  sig- 
nificance level-by-level  than  those  of  either  COEF-1  or 
C0EF-2S  coefficient  sets  when  applied  to  the  same  data. 
These  results,  as  much  as  any  individual  ship  results,  indi- 
cated the  need  to  consider,  initially,  a  lat itudinally- 
stratified  set  of  predictors,  with  the  boundary  set  at  <j)_>50N. 
In  addition,  COEF-1  regressions  consider  the  alternate  hy- 
pothesis, namely,  the  case  of  no  latitudinal  stratification. 
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Note  that  like  results  for  the  remaining  nine  ships 
of  the  OSV-2  sample  are  listed  in  Appendix  B  in  the  order 
specified  by  Table  3. 

2.   For  the  Determination  of  a  "Best"  Set  of  Coefficients 
a.   Level-by-Level 

Water  vapor  estimates  were  made  by  COEF-1  and 
C0EF-2N,S  for  the  entire  independent  OSV-2  sample.   The  re- 
sults are  to  be  tested  level-by-level  as  follows: 

i.   COEF-1  vs.  C0EF-2N  for  latitudes  cJ)>_50N 
ii.   COEF-1  vs.  C0EF-2S  for  latitudes  <})<50N 
iii.   COEF-1  vs.  COEF-2   for  all  <j> 
The  statistical  results  for  the  three  comparisons  are  listed 
in  Tables  6  through  8,  respectively.   As  in  Section  IV. A  the 

F  value  for  a  1%  significance  level  was  taken  for  the  most 
c 

stringent  case,  i.e.  case  (ii)  was 
F  (1,250)  =  6.75  . 

The  summary  of  F-statistics  in  Table  7  for  <J)_>50N 
shows  strong  predictive  capability  at  all  levels  in  both 
COEF-1  and  C0EF-2N,  with  the  latter  demonstrating  slightly 
higher  verification  (F-value)  statistics.   The  results  of 
C0EF-2S  were  not  as  strong  when  compared  to  COEF-1.   Table  8 
shows  slightly  better  verifications  with  COEF-1  from  level  1 
through  11,  then  a  reversal  for  the  remaining  eight  levels. 
Note  that  at  level  9,  for  the  C0EF-2S  tabulation,  the  value 
R   =0,  but  also  F<0.   The  Rp^-value,  as  defined  by  Eq.  (22) 
takes  on  negative  values  whenever  the  mean  square  error  of 
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the  estimation  scheme  exceeds  the  sample  variance.   For  this 
case,  the  sample-mean  must  be  taken  as  the  best  prediction 
scheme,  and  R-p-p-p  is  reset  to  0  in  obtaining  a  realistic 
value  of  RFp-p  •   However  the  F-value  was  left  as  calculated 
by  Eq.  (23),  with  k=l  and  N-l  degrees  of  freedom.   The  sta- 
tistical comparison  of  COEF-1  and  COEF-2,  listed  in  Table  9, 
for  the  entire  OSV-2  independent  sample  was  a  composite  of 
the  results  of  Tables  7  and  8.   Both  estimation  methods 
demonstrated  strong  prediction  qualities,  as  evidenced  by 
the  large  F-values.   As  in  the  regression  analyses  in  Section 
IV. A,  the  comparison  of  the  performance  of  COEF-1  and  COEF-2, 
Table  9,  is  quite  close  over  most  of  the  19  levels.   Hence, 
there  was  no  conclusive  evidence  supporting  a  decision  to 
select  a  ,:best"  set  of  coefficients  from  either  COEF-1  or 
COEF-2.   This  decision  was  deferred  to  the  outcome  of  the 
regression  results  relating  observed  to  estimated  values  of 
precipitable  water  vapor. 

b.   Precipitable  Water 

Precipitable  water  calculations  (U)  were  made  on 
the  OSV-2  sample  soundings  according  to  Eq.  (7).   Similar 
application  of  (7)  to  the  19-level  OSV-2  mixing  ratio  speci- 
fications generated  by  COEF-1  and  the  composite  set  of 
COEF-2  equations  produced  two  modes  of  estimated  U-values 
for  each  of  the  829  soundings,  depending  upon  the  latitude 
of  the  sounding. 

The  simple  linear  regression  analyses  were  per- 
formed on  each  set  of  U,  U  using  BMD03R.   Each  regression, 


42 


CN 
I 

En 

w 
o 


CO 

o 

e 
o 
o 

c 

CO 


I 

W 

O 

c_> 

>s 
x> 

c 
o 

•H 
4-1 
CO 

a 


co 


c 

C 

ex 

0) 

c 

•H 


c 
o 
w 

•H 
>-i 
CO 
Cu 
E 
O 

o 


CO  • 

o  -e- 

•H 

•U  iH 

W  i-H 

•H  CO 
4-1 

CO  U 

4J  O 

to  iw 


on 

4) 
iH 

CO 
H 


in 

o 

m 

O 

p- 

f—< 

in 

CO 

r- 

C\I 

o 

>*• 

P- 

o 

oo 

■— i 

sO 

(NJ 

sO 

«o 

—4 

sO 

— i 

o 

in 

in 

CNI 

•—4 

r»- 

vt 

ro 

^O 

4- 

•* 

in 

o 

in 

1—4 

>o 

O 

>t 

CNJ 

o 

h- 

*t 

r- 

in 

CNJ 

CNJ 

m 

t-i 

CO 

00 

CNI 

o 

ro 

oo 

u_ 

«t 

CO 

CNJ 

CO 

—1 

O 

r— 4 

^t 

t-4 

<f 

0> 

CNJ 

o 

CNJ 

-o 

O 

p- 

on 

^-4 

o 

r-1 

— 

— - 

CNJ 

o 

st 

CO 

P- 

CNJ 

CO 

o 

<t 

1—4 

in 

CO 

CO 

on 

ON 

CN 
I 

LO 

o 

O 

o 

o 

vO 

r»- 

» 

r-4 

O 

in 

m 

in 

vO 

CO 

ro 

—4 

CO 
CNJ 

p- 
in 

ON 
sO 

En 

Ed 

O 

U 

u. 

CO 

■n 

ro 

c- 4 

ON 

o 

m 

co 

in 

ON 

r-4 

co 

vO 

«d- 

<st 

CO 

ON 

m 

s0 

Ed 

u_ 

vO 

CM 

CNJ 

st 

st 

st 

r- 

h- 

^i- 

in 

in 

st- 

o 

CNJ 

CO 

o 

ON 

in 

in 

H 

LU 

st 

in 

in 

in 

m 

r- 

00 

O 

^-4 

in 

sf 

—4 

("0 

in 

i—i 

ON 

p- 

CO 

st 

M 

on 

Ct 

o 

sO 

sO 

s0 

so 

sG 

o 

s0 

<f 

vO 

>J3 

n0 

m0 

s0 

r»- 

p- 

00 

a* 

O' 

O 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

u 

ON 

~ 

p- 

st 

ro 

o 

o 

ON 

ro 

s0 

r^ 

o 

CO 

^ 

o 

ON 

sO 

p- 

m 

•    • 

o 

tn 

CO 

CNJ 

CNJ 

o 

rO 

m 

CO 

o 

r- 

CO 

•4- 

■st- 

•st 

CO 

CO 

co 

p- 

QQC 

HOC 

— . 

—4 

CNJ 

ro 

st 

>o 

CO 

-t 

f— 4 

«* 

CO 

in 

p- 

—* 

ON 

CO 

o 

O 

sO 

00  LU 

»-4 

ro 

ro 

CO 

m 

-o 

CO 

CO 

CO 

CO 

p- 

p- 

in 

o 

CO 

CNJ 

O 

CNI 

in 

ON 

f-4 

vO 

in 

m 

*t 

l»»- 

•—4 

in 

p- 

ro 

rr, 

p- 

f\J 

co 

1—4 

in 

ro 

ro 

vO 

■4" 

o 

CO 

•—4 

p- 

00 

i—4 

CO 

o 

sO 

ro 

>o 

o 

r— 4 

in 

—4 

p- 

•—4 

in 

in 

■—4 

o 

^_ 

n£) 

ro 

ON 

sO 

CNJ 

O 

s0 

LL 

ro 

CO 

a> 

ro 

<0 

i— i 

o 

^ 

o 

CNJ 

p- 

r* 

P- 

CNJ 

st 

r- 

CNJ 

sO 

CNI 

n0 

r- 

f»- 

CO 

CO 

ON 

in 

r- 

•—i 

CNI 

CO 

r- 

i—4 

00 

0> 

p- 

ro 

CNJ 

CNI 

in 

in 

in 

in 

in 

o 

h- 

h- 

p- 

v0 

in 

«* 

in 

in 

P- 

r\j 

•—4 

r— 
CNJ 

in 
in 

CNJ 

P- 

t-l 

u. 

vO 

«t 

co 

CNJ 

o 

O 

vO 

<t- 

00 

(NJ 

st 

ON 

in 

oo 

r-4 

CNJ 

—4 

sO 

^ 

1 

Pu 

Ed 
O 

LL 

•o 

— < 

•—4 

ro 

st 

st 

on 

in 

0* 

in 

■4" 

<h 

o 

CNJ 

o 

ON 

sO 

C\.l 

P- 

LU 

rO 

<r 

*t 

<t 

<* 

P- 

CO 

a> 

r- 

in 

•d- 

O 

CNJ 

•J- 

o 

P~ 

P- 

ro 

•t 

U. 

•O 

v0 

-o 

O 

o 

s0 

o 

v0 

vO 

o 

O 

>o 

>0 

nQ 

P- 

r>- 

co 

ON 

ON 

C_! 

O 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O 

o 

O 

<"0 

o 

CO 

o 

o 

CO 

>d- 

in 

CO 

o 

o 

ON 

CNJ 

ON 

r>t 

ro 

co 

st 

•    • 

— < 

in 

«t 

CNJ 

CNJ 

o 

CO 

o 

o 

o 

CO 

sj- 

*— 1 

CO 

o 

o 

in 

CO 

s0 

Co: 

r-rX 

T-t 

1—4 

CNJ 

co 

■tf 

o 

CO 

<)■ 

in 

vf 

CO 

•0 

CO 

CNJ 

—4 

i—4 

r-4 

CNJ 

in 

OOLU 

i—4 

CNJ 

CO 

CO 

m 

o 

CO 

ON 

ON 

CO 

^ 

p- 

«\ 

o 

ro 

CO 

o 

rr\ 

i-4 

P- 

CO 

CNJ 

*d- 

ro 

in 

-o 

j   1 

p- 

r-4 

in 

^-4 

o 

•     t 

>t 

o 

•—i 

<\j 

o 

—t 

in 

CO 

CT 

1—4 

r- 

00 

o 

<t 

m 

CNJ 

i—i 

sO 

o 

Q> 

1— LU 

1—1 

rsj 

ro 

^t 

in 

CO 

•—4 

o 

Nf 

in 

O 

o 

•£> 

p- 

p- 

in 

O 

o 

sO 

00  CJ 

^H 

CNJ 

ro 

vl- 

in 

r- 

CO 

o 

CNJ 

st 

sQ 

o 

ro 

o 

•—4 

i-4 

f— 4 

—4 

CNJ 

Psl 

CNJ 

CO 

II 

2 

o 

CNI 

p- 

IT. 

m 

CO 

n0 

vO 

in 

m 

o 

o 

p- 

•—4 

rO 

CNJ 

st 

r»- 

CNJ 

<: 

cr> 

f- 

CVi 

CO 

-0 

tn 

r— 1 

r*- 

on 

sC 

sj- 

o 

ro 

sO 

O 

rs- 

r~ 1 

CNJ 

p- 

LU 

LU 

rvj 

5: 

.— 1 

f\J 

<*• 

in 

p- 

o 

J" 

CNJ 

r» 

o 

0N 

in 

— 

ro 

—4 

s± 

ro 

p- 

ro 

*-* 

i— 4 

t— 4 

CNJ 

ro 

in 

in 

CO 

i—4 

>cf 

ON 

m 

st 

p~> 

o- 

oO 

^-4 

r-4 

i—4 

CNJ 

ro 

st 

m 

U_ 

_J 

_i 

CL 

LU 

s 

> 

.—1 

cnj 

ro 

st 

in 

>o 

r- 

co 

O 

O 

—- ( 

CNJ 

CO 

<1- 

in 

s0 

P- 

CO 

ON 

< 

LU 

oo 

-J 

43 


identical  to  the  MLR  with  only  one  independent  variable,  was 
performed  with  U  as  the  dependent  variable  and  the  appro- 
priate  U  as  the  independent  variable.   These  regressions 
provided  significance  measurements  of  the  19-level  predic- 
tions by  a  single  vertically-scaled  F-statistic  with  1  and 
N-2  degrees  of  freedom. 

The  statistical  results  from  the  precipitable 
water  vapor  analyses  are  listed  in  Table  10  and  provided  a 
final  comparison  between  the  COEF-1  and  COEF-2  specification 
systems.   The  critical  F-value, 

F   =  7.88 
c 

corresponding  to  the  sample  of  N=251  soundings  (4><50N)  pro- 
vided the  most  stringent  F-statistic  significance  require- 
ments namely  at  the  1%  level.   As  seen  in  Table  10,  all 
calculated  F-values  greatly  exceeded  the  critical  value, 
thus  indicating,  to  at  least  the  1%  level  of  significance, 
that  none  of  the  specifications  of  precipitable  water  could 
have  occurred  by  chance. 

The  precipitable  water  results  were  analogous  to 
those  in  Tables  7  through  9,  and  as  seen  in  Table  10,  the 
single  COEF-1  set  of  equations  yielded  the  most  effective 
specifications.   Considering  the  conditions  set  down  in  Sec- 
tion IV.B.l  requiring  considerably  better  specification  by 
COEF-2  over  COEF-1,  and  with  the  results  of  Sections  IV. A 
and  IV. B,  the  decision  was  made  at  this  point  to  eliminate 
from  further  consideration  the  model  of  a  lat itudinally 
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differentiated   set   of   coefficients.      Thus   the   remainder   of 
this   research    incorporated   only   the   composite  COEF-1   set   of 
equations    for   comparison  with   the  WC   specification  method. 

Table   10.      Precipitable  water  vapor   regression  statistics   of  U  on  U  by 
COEF-1   and  by   the  appropriate  COEF-2   scheme  applied    to   inde- 
pendent  OSV-2   sample. 


COEF-1 

C0EF-2(N,S,  Composite) 

AREA 

MEAN 

STD 
DEV 

STD 
ERR 

R 

F 

STD 
ERR 

<t»50 

.998 

.408 

.2344 

.8191 

1170.1 

.2418   .8054  1063.8 

(N=578) 

4><50 

1.73 

.674 

.3847 

.8218 

517.9 

.3902   .8161   496.4 

(N=251) 

All  <J> 

1.22 

.606 

.2881 

.8800 

2838.6 

.2947   .8741  2677.2 

(N=829) 



<■    ■  ■' 

Operational    implementation   of    a   single   set    of 
coefficients   over   a  more   complex   set    eliminated   the   require- 
ment   of    additional   computer   core   allocation   for    storage   of    a 
second   set    of   coefficients   and   the   necessary    logic   state- 
ments  to   enact    the   stratification   differentiation.      Clearly, 
the   current    results   do   not   warrant   use   of   the  more   complex 
COEF-2   prediction    scheme. 

3.       Comparison   With   the   WC   Predictors 
a.      Level-by-Level 

In  this  section  comparison  of  estimates  produced 
by  the  "best"  COEF-1  coefficient  set  was  made  with  those  es- 
timates   from   the  WC   predictions.      Comparison  was   conducted 
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in  accordance  with  the  precedures  discussed  in  Section 

IV.B.l.a.   As  in  previous  comparisons,  the  most  stringent 

conditions  were  adhered  to  for  determination  of  an  F  -value, 

c       ' 

in  this  case 

F  (1,250)   =  6.75  . 

Statistical  comparison  of  the  two  specification  methods  are 
listed  in  Tables  11,  12,  and  13  for  the  latitudinally  strati- 
fied subsamples  (0SV-2N,  0SV-2S)  and  for  the  total  OSV-2 
sample,  respectively. 

Without  exception,  the  COEF-1  specifications 
were  statistically  superior  to  those  produced  by  the  WC  co- 
efficients at  all  19  J-levels  in  each  of  the  three  tables. 
The  WC  specifications  do  demonstrate  considerable  predicta- 
bility in  the  lower  levels  of  the  atmosphere,  however,  the 
F-values  indicate  poor  results  in  the  highest  two  or  three 
levels.   Some  possible  explanations  will  be  offered  in  Sec- 
tion V  for  the  breakdown  of  the  upper  level  specifications 
by  the  WC  equations.   It  is  clear  from  the  results  listed  in 
Tables  11  through  13  that  the  COEF-1  specifications  are  su- 
perior to  the  WC  predictions  for  the  level-by-level  analysis 
b.   Precipitable  Water 

A  final  comparison  of  the  WC  and  COEF-1  specifi- 
cation methods  was  made  using  precipitable  water  vapor  cal- 
culations as  in  Section  IV.B.l.b.   Table  14  lists  the 
statistics  resulting  from  the  vertically-scaled  evaluation 
of  the  two  water  vapor  estimates  for  each  of  the  OSV-2 
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stratifications.   The  statistics  were  generated  from  two 
simple  linear  regressions  by  BMD03R  of  U  on  each  of  the  U- 
estimates.   The  last  line  of  Table  14  indicates  a  definite 
superiority  of  the  COEF-1  equations  over  the  WC  equations 
for  the  independent  OSV-2  data  sample. 

Table  14.   Precipitable  water  vapor  results  by  COEF-1  and  modified  WC 
coefficients  applied  to  independent  OSV-2  sample. 


COEF-1 

WC 

AREA 

MEAN 

STD 
DEV 

STD 

ERR     R 

F 

STD 
ERR 

R 

F 

cp>50 

.998 

.408 

.2344   .8191 

1170.1 

.2658 

.7587 

781.1 

(N=578) 

cp<50 

1.73 

.674 

.3847   .8218 

517.9 

.5324 

.7297 

283.5 

(N=251) 

All  <|> 

1.22 

.606 

.2881   .8800 

2838.6 

.3393 

.8289 

1816.1 

(N=829) 

C.   ISLAND  AND  COASTAL  (NESS)  DATA 
1.   Level-by- Level 

The  160  NESS  island  and  coastal  March  soundings  pro- 
vided a  test  sample  that  focussed,  to  a  degree,  on  evaluation 
of  the  possible  significance  of  winter  effects  on  W-profile. 
Table  15  presents  the  statistics  obtained  by  the  first  ef- 
forts at  estimating  the  W-profiles  by  COEF-1  and  WC  coeffi- 
cients.  In  making  these  first  W-specif ications  all  160 
soundings  (with  observed  W-values  up  to  at  least  510  mb) 
were  accepted  as  valid  data.   The  results  presented  in  Table 
15  were  at  best  discouraging  for  both  sets  of  coefficients. 
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The  negative  F-values,  starting  at  level  13  for  the  COEF-1 
specifications  and  at  level  9  for  those  by  the  WC  method, 
indicated  considerably  greater  effectiveness  would  be  achieved 
by  making  predictions  with  the  sample  mean  values  of  the  mix- 
ing ratios. 

It  was  noted  that  a  number  of  the  NESS  soundings 
should  not  have  been  considered  reliable  after  reviewing  the 
actual  W-profiles.   Unreliable  soundings  were  identified 
when  the  reported  moisture  profiles  indicated  RH<20%  in  the 
lower  levels  of  the  soundings  and  persisted  with  increasing 
ascent.   It  was  felt  that  the  low  moisture  reports  were  due 
mostly  to  a  failure  of  the  moisture  element  on  the  radio- 
sonde instrument. 

If  a  sounding  had  RH  considerably  higher  than  20%  at 
all  levels,  its  moisture  profile  was  considered  reliable  and 
the  sounding  was  retained.   In  some  cases  the  relative  hu- 
midity decreased  to  values  lower  than  20%  in  the  mid-tropo- 
sphere but  then  recovered  to  exceed  20%  at  higher  levels . 
Such  soundings  were  also  retained  in  a  revised  test  set.   If 
the  RH  persisted  at  less  than  20%  throughout  most  of  the 
troposphere,  including  the  upper  regions  where  motor-boating 
is  not  infrequent,  the  sounding  was  considered  unreliable 
and  discarded  (see  Table  16). 
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Table  16.  Examples  of  three  types  of  soundings  categorized 
by  RH  profiles  extending  up  to  level  10  and  con- 
sidered for  possible  inclusion  in  the  NESS  sample 


evel 

Reliable 

(%  RH) 
Recovered 

10 

72 

71 

11 

67 

63 

12 

62 

23 

13 

45 

11 

14 

42 

9 

15 

43 

11 

16 

42 

16 

17 

37 

31 

18 

57 

29 

19 

79 

38 

Unreliable 


15 
16 
16 
17 
18 
21 
20 
18 
67 
76 


After  screening  the  NESS  soundings  using  the  method 
of  Table  16,  a  total  of  124  soundings  remained  for  testing. 
Table  17  lists  the  statistical  results  of  the  level-by-level 
W-specif ications  by  the  COEF-1  and  WC  equation  sets.   Since 
the  number  of  soundings  with  observed  W-values  drop  off 
sharply  above  level  10,  the  critical  F-value  was  arbitrarily 
taken  for  N=70  occurring  at  level  7  and 


F   =7.01 
c 


This  value  of  F   is  considerably  more  stringent  than  neces- 
sary  for  specifications  at  levels  below  level  J=7,  but  as 
Table  17  demonstrates,  this  is  a  reasonable  restriction. 
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Examination  of  the  F-values  of  Table  17  indicated 
that  the  WC  coefficients  performed  more  effectively  in  the 
lower  levels  of  the  atmosphere.   Above  level  16,  the  COEF-1 
set  of  equations  produced  statistically  superior  estimates 
up  to  level  4,  where  both  estimation  procedures  produce  neg- 
ative R-^-p-p  and  F-values.   The  small  number  of  soundings  con- 
taining data  in  the  last  four  levels  was  considered  responsible 
for  these  unusual  statistical  outcomes  and  it  was  impossible 
to  draw  any  useful  conclusion  regarding  the  comparative 
upper-level  specificity  by  the  two  estimation  systems. 
2.   Precipitable  Water  Vapor 

Precipitable  water  vapor  integrations  were  performed 
using  Eq.  (7)  on  each  of  the  124  NESS  soundings,  incorpo- 
rating all  available  J-levels  for  which  mixing  ratio  values 
existed,  J<_9.   Additionally,  U-values  were  calculated  using 
only  the  corresponding  levels  for  the  associated  U-integra- 
tions  up  to  level  J. 

The  statistical  results,  as  produced  by  a  BMD03R 
simple  linear  regression  of  U  on  each  of  the  two  U's,  are 
displayed  in  Table  18. 

Table  18.   Statistical  results  from  BMD03R  simple  linear 

regressions  on  precipitable  water  vapor  calcula- 
tions obtained  from  the  124  NESS  observed  and 
COEF-1,  WC-estimated  mixing  ratio-values. 


COEF-1 

WC 

MEAN 

STD 
DEV 

STD 
ERR 

R 

F 

STD 
ERR 

R 

F 

1.44 

1.1 

.4292 

.9218 

690.1 

.4850 

.8990 

514.0 
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The  critical  F-value 
F  (1,122)  =  6.84 

was  compared  to  the  calculated  F-values  of  Table  18  to  eval- 
uate the  significance  of  the  vertically  scaled  mixing  ratio 
specifications.   It  was  readily  apparent  that  the  U-estimates 
both  by  the  COEF-1  and  WC  systems  demonstrated  high  statis- 
tical significance  and  that  the  COEF-1  estimates  indicated 
somewhat  higher  predictability. 

Reviewing  the  results  of  Table  18  along  with  those  of 
Table  17,  it  became  obvious  that  the  near-surface  mixing 
ratio  values  dominate  the  precipitable  water  vapor  integra- 
tion, and  that  the  upper  level  breakdown  of  the  WC  specifi- 
cations, as  evidenced  by  Table  17,  was  masked  by  this 
dominance.   Note  that  the  W-values  of  the  dependent  and  in- 
dependent OSV  data  sets  were  determined  by  Smith's  power  law, 
Eq.  (6),  for  levels  above  J=9.   On  the  other  hand,  the  W- 
values  for  the  NESS  data  set  were  taken  solely  as  indicated 
by  the  radiosonde  instrument,  as  were  the  data  used  by  WC  in 
determining  their  regression  coefficients.   Application  of 
the  power  law  to  the  OSV  data  tended  to  smooth  the  W-profile 
at  higher  levels  and  yielded  a  degree  of  consistency  in  the 
OSV  W-profiles  in  the  highest  5  prediction  levels  (P£300  mb ) , 
giving  small  W-values  over  the  OSV's  at  levels  higher  than 
J=5  as  compared  with  those  resulting  from  the  NESS  sounding 
samples  (Table  17).   It  is  felt  that  the  high  values  of  W  in 
the  NESS  data  could  be  attributed  largely  to  inaccuracies  in 
the  radiosonde  observational  data  at  the  higher  levels. 
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V.   CONCLUSIONS 

This  thesis  has  reinforced  the  findings  of  Weinreb  and 
Crosby  (WC)  in  providing  supporting  evidence  that  mixing 
ratios  can  be  described  by  a  linear  combination  of  a  desig- 
nated 11-element  temperature  T(P)-prof ile ,  which  is  equiva- 
lent to  a  profile  of  W  (P).   Further,  it  has  been  shown  that 
any  dependence  of  the  relationship  between  temperature  and 
mixing  ratio  on  latitude  is  insignificant,  and  that  one 
single  set  of  specification  equations  can  be  employed  over 
all  latitudes. 

It  was  evident  from  the  statistical  results  presented  in 
Section  IV. B. 2  that  the  performance  of  the  coefficient  set 
(COEF-1)  was  superior  to  that  of  the  WC  coefficient  set  over 
the  Ocean  Station  Vessels.   However,  the  causes  of  the  im- 
proved specifications  were  not  regarded  as  a  result  of  a 
land-mass  (coastal)  or  seasonal  effect.   Rather,  inaccuracies 
inherent  in  radiosonde  moisture  data  in  the  higher  level 
atmospheric  structure  may  have  affected  the  definition  of 
the  WC  coefficients,  which  were  derived  by  moisture  data  of 
accuracy  comparable  to  that  of  the  NESS  data.   This  reason- 
ing is  supported  by  the  independent  tests  of  the  COEF-1  and 
WC  specifications  on  the  NESS  data. 

Efforts  to  produce  an  accurate  prediction  scheme  by  MLR 
techniques,  as  demonstrated  by  this  work,  appear  extremely 
promising.   Effective  prediction  procedures  seem  not  to  be 
associated  with  seasonal,  latitudinal  or  land-mass  (coastal) 
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effects.   Rather,  effective  prediction  appears  to  depend 
upon  an  accurate  data  base  of  sufficiently  good  quality  from 
which  to  define  accurate  coefficient  sets  over  the  depth 
1000-200  mb.   The  results  of  this  research  indicate  that  the 
OSV  data  provide  such  a  data  base,  and  that  mixing  ratio  es- 
timates by  a  set  of  equations  of  the  form  derived  here  pro- 
vide a  representative  profile  of  W(P)  up  to  approximately 
200  mb.   The  high  level  of  predictability  of  the  precipitable 
water  vapor  resulting  from  the  vertical  integration  of  mixing 
ratio  estimates  provided  substantial  evidence  supporting  the 
conclusion  that  the  OSV  data  define  highly  effective  regres- 
sion equations. 


58 


APPENDIX  A 

This  Appendix  contains  a  listing  of  the  three  sets  of 
regression  coefficients  COEF-1,  C0EF-2N,  and  C0EF-2S  deter- 
mined by  B.MD03R  from  the  OSV-1  ship  sample.   A  fourth  list- 
ing contains  the  selection  of  coefficients  extracted  from 
the  WC  set  described  in  Section  I. 

Each  of  the  four  sets  of  coefficients  was  applied  accord- 
ing to  the  following  matrix  operation  to  produce  the  appro- 
priate mixing  ratio  estimates  at  level  J: 
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APPENDIX  B 

Appendix  B  lists  the  level-by-level  statistical  results 
obtained  by  application  of  coefficient  sets  COEF-1,  C0EF-2N, 
and  C0EF-2S  to  the  individual  OSV-2  soundings.   The  listings 
are  arranged  in  the  order  specified  by  Table  3,  only  that 
Ship  Station  P  is  not  included  here  (see  Table  6). 
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